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ABSTRACT 

We examine isothermal dark matter halos in hydrostatic equilibrium with a"A- 
field", or cosmological constant A = Q\p CIit c 2 , where Ha — 0.7, and p crit is the 
present value of the critical density with h ~ 0.65. Modeling cold dark matter as a 
self-gravitating Maxwell-Boltzmann gas, the Newtonian limit of General Relativity 
yields equilibrium equations that are different from those arising by merely coupling 
an "isothermal sphere" to the A-field within a Newtonian framework. Using the condi- 
tions for the existence and stability of circular geodesic orbits, the numerical solutions 
of the equilibrium equations (Newtonian and Newtonian limit) show the existence of 
(I) an "isothermal region" (0 < r < r%), where circular orbits are stable and all vari- 
ables behave almost identically to those of an isothermal sphere; (II) an "asymptotic 
region" (r > r±) dominated by the A-field, where the Newtonian potential oscillates 
and circular orbits only exist in disconnected patches of the domain of r; (III) a "tran- 
sition region" (r% < r < r±) between (I) and (II), where circular orbits exist but are 
unstable. We also find that no stable configuration can exist with central density, p c , 
smaller than 2 A, hence any galactic haloes which virialized at z < 30 in a A-CDM 
cosmological background must have central densities of p c > 0.008 M Q /pc 3 , in inter- 
esting agreement with rotation curve studies of dwarf galaxies. Since marks the 
largest radius of a stable circular orbit, it provides a characteristic boundary or "cut 
off" maximal radius for isothermal spheres in equilibrium with a A-field. For current 
estimates of p c and velocity dispersion of virialized galactic structures, this cut off 
scale ranges from 90 kpc for dwarf galaxies, up to to 3 Mpc for large galaxies and 22 
Mpc for clusters. In a purely Newtonian framework these length scales are about 10% 
smaller, though in either case ri is between five and seven times larger than physical 
cut off scales of isothermal halos, such as the virialization radius or the critical radius 
for the onset of Antonov instability. These results indicate that the effects of the A- 
field can be safely ignored in studies of virialized structures, but could be significant 
in the study of structure formation models and the dynamics of superclusters still in 
the linear regime or of gravitational clustering at large scales (r w 30 Mpc). 

Key words: gravitation - relativity - galaxies: haloes - galaxies: formation - cos- 
mology: theory 



1 INTRODUCTION 

Recent observational data from type la supernovae and 
CMB anisotropics suggests that large scale cosmic dynamics 
are dominated by a repulsive type of matter-energy density, 
generically known as a "dark energy". In its most simple 
empirical formulation, dark energy can be described as a 
"A-field" associated with the "cosmological constant" . More 
sophisticated models of dark energy take the form of scalar 



field sources in which the A-field corresponds to the limit 
of vanishing scalar field gradients. This paper addresses an 
important issue in this context, namely: to examine if this 
type of repulsive interaction has any influence on the equi- 
librium of virialized galactic halo structures. Various aspects 
regarding this question have been considered previously in 
the literature. Lahav et al. (1991) have argued, using Newto- 
nian spherical infall models (or "top hat" models), that the 
effects of a cosmological constant can be neglected in the 
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study of structure formation. Further improvements of "top 
hat" models incorporating A and/or quintessence-like fields 
have been considered by Wang & Steinhardt (1998), Iliev & 
Shapiro (2001) and Lokas & Hoffman (2001). Constraints on 
A at the level of clusters and superclusters have been placed 
by means of optical gravitational lens surveys (Quast & Hel- 
big 1999) and peculiar velocities of galaxies (Zehavi & Dekel 
1999). Other papers (e.g. Axenides et al. 2000, Nowakowski 
2001 and Nowakowski et al. 2001) have studied the effect 
of A on galactic structures within the framework of Gen- 
eral Relativity, though looking at this effect not in the self- 
gravitating sources themselves, but through the dynamics of 
test observers in a Schwarzschild-deSitter field that is sup- 
posed to be the "external" field of these sources. 

In this paper we complement and extend previous 
work by studying a spherically symmetric, self-gravitating, 
Maxwell-Boltzmann (MB) gas in hydrostatic equilibrium 
with a A-field given as a "cosmological constant" , subject to 
observational constraints. Although galactic halo structures 
are characterized by Newtonian energies and velocities, we 
find it to be more correct to examine these sources within 
General Relativity and then obtain the adequate ("weak 
field") Newtonian limit. The reason for this approach follows 
from the fact that in Newtonian physics the mass density is 
the only source of gravitation and so a non-gravitational 
interaction, such as the A-field, does not contribute to the 
"effective gravitational mass". In General Relativity, on the 
other hand, all forms of matter-energy count as sources of 
gravitation (curvature of spacetime). Therefore, even in the 
Newtonian limit, the matter-energy contribution from the 
A-field could be comparable (or at least might not be negli- 
gible in comparison) to rest-mass density of the MB gas, so 
that the hydrostatic equations furnished by the Newtonian 
limit might not coincide with those "naively" obtained from 
a Newtonian framework. As we show in this paper, for a MB 
gas and a constant A-field, the Newtonian limit yields a dif- 
ferent set of equilibrium equations from those of Newtonian 
theory. However, we cannot know a priori if this difference 
leads to significant or negligible effects until we solve the 
equilibrium equations and compare their solutions in both 
cases. Such ambiguity regarding the Newtonian limit would 
not arise if we were only dealing with a MB gas (a thermal 
source), this is so because in this case the internal energy 
(the only contribution to matter-energy besides rest-mass) 
is negligible in comparison with rest-mass, leading to the 
"isothermal sphere" as the unequivocal Newtonian limit. 

It is a well known result that within the hierarchical 
model of structure formation (White & Rees 1978) the final 
halo structure which results from the repeated merger events 
building up a galaxy is well represented not by an isother- 
mal structure, but by a centrally divergent steep inner halo 
profile, first noted by Navarro, Frenk & White (1997), and 
since confirmed by numerous authors. However, controversy 
still surrounds this result in as much as direct dynamical 
studies of galactic rotation curves, across wide dynamical 
ranges in galaxy size, appear to indicate the presence of a 
constant density core in the dark haloes of galaxies (e.g. 
Burkert & Silk 1997, de Blok, McGaugh & Rubin 2001, Fir- 
mani et al. 2001, among others). This being the case, we 
think it worthwhile to study in depth the modifications a 
simple isothermal halo suffers in the presence of a domi- 
nant A component, even though this results may prove to 



be merely indicative of a case where not just a stable halo 
is being treated, but including the full merger-driven for- 
mation of such a structure. Whatever the correct structure 
scenario eventually turns out to be, a detailed understand- 
ing of the behaviour of the quintessential galactic halo toy 
model -the isothermal halo- in the presence of an important 
cosmological A term, will surely prove valuable. 

The organization of the paper is described below. Sec- 
tions II, III and IV respectively provide the general relativis- 
tic formulation for a non-relativistic MB gas, a A-field as 
a zero gradient limit of a quintessence-like scalar field and 
the equations of hydrostatic equilibrium for these sources. 
We examine the Newtonian limit of the models in section V, 
showing that in the equilibrium equations in this limit do 
not coincide with the corresponding Newtonian equations 
for the same sources. In section VI we discuss the cut off 
length scales, n, r2, introduced by the A-field, comparing 
them in section VII (by means of rough qualitative argu- 
ments) with physically motivated length scales of virialized 
halo structures (virialization radius and critical radius asso- 
ciated with Antonov instability). In section VIII we review 
briefly the equilibrium equations for the King halos, while 
in section IX we use the numerical integration of the vari- 
ous sets of equilibrium equations in order to verify the issues 
discussed analytically and qualitatively in previous sections: 
the difference between the Newtonian limit of General Rela- 
tivity and a pure Newtonian framework and the fact that ri 
is 5-7 times larger than physically motivated length scales. 
We briefly summarize these results in section X. The paper 
contains 9 figures obtained from the numerical solution of 
the equilibrium equations. 



2 THE NON-RELATIVISTIC 

MAXWELL-BOLTZMANN GAS 

A relativistic MB gas follows from a Kinetic Theory for- 
malism associated with the equilibrium MB distribution (de 
Groot et al. 1980). For gas particles of mass m, the param- 
eter of "relativity coldness" is defined as 



k B T> 



(1) 



where k B is Boltzmann constant and T is the temperature. 
The non-relativistic MB gas corresponds to the limit j3 ^> 1, 
leading to the following forms for the energy density s and 
pressure p 



n ^ 



2 . 3 U 

inc + - k E 



P 



nk„ T = 







with the particle number density, n, given by 

3 

exp I 



n = 



mc 



(2) 
(3) 

(4) 



where 1-kTi is Planck's constant and ^ is the Gibbs function 
per particle (de Groot et al. 1980), the chemical potential. 
The MB distribution is a solution of the Einstein-Boltzmann 
equations, irrespective of whether the gas is collisional or col- 
lisionless, its applicability is only restricted by the criterion 
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of "non-degeneracy" or "dilute occupancy" (Pathria 1972): 
exp(fif3/mc 2 ) <JC 1. For typical galactic halo variables, this 
criterion holds for particle masses complying with m >• 60 
eV (Cabral-Rossetti et al. 2002). 

Under General Relativity, a spacetime whose source is a 
self gravitating MB gas can be described by the Momentum- 
Energy tensor of a perfect fluid 

T,ab a b . 7 ab /p-\ 

MB = e u u + ph , (5) 

where e, p follows from equations (2)-(4) and h ab = u a u b + 
g ab . Such a spacetime must admit a timelike Killing vector 
field (Maartens 1996) 



p a =(3u a , 
so that 

Ua-h b a {\np), b = 0, 



(6) 



(7) 

where ii a = u :b u b is the 4-acceleration. Condition (6) implies 
that the spacetime must be stationary (or static if vorticity 
vanishes), while condition (7) is the "Tolman Law" indi- 
cating that equilibrium of a MB gas in General Relativity 
requires the existence of a nonzero temperature gradient, a 
specific spacelike gradient that balances the 4-acceleration 
associated with a stationary frame. In a comoving represen- 
tation, condition (7) becomes T^J—gn — const. 



3 QUINTESSENCE AND THE A FIELD 

"Dark energy" matter-energy sources are usually rep- 
resented as a quintessence-like scalar field <f> with a 
Momentum-Energy tensor 

1 , 



T<ab 
O — 



■ + V(4>) 



(8) 



where V(<j)) is the scalar field potential. Since we are not 
interested in a cosmological quintessence source, but rather 
in its gravitational interaction with a MB gas at the galactic 
and galactic cluster scale, it is reasonable to assume that 
the gradients 0, a are negligible, hence we can approximate 
equation (8) by a "cosmological constant" or a A-field given 
by 

T<f « -V(<t> )g ab = -Ag ab , (9) 

equivalent to a source with matter-energy density 

A = SI a Pc?it c 2 , 
p crit = 1.88 x 10~ 29 h 2 gm/cm 3 (10) 

where, according to observational evidence, we can take 
Qa — 0.7 and h ~ 0.65. If there is only gravitational inter- 
action between the MB gas and the A field, the conditions 



Q ;f> 



0, 



rpab 
1 MB ; b 



0, 



(11) 



must hold separately. This condition is trivially satisfied for 
Tq given by equation (9). 



4 THE FIELD EQUATIONS 

In order to model an equilibrium galactic halo structure we 
will consider spherically symmetric static spacetimes, com- 



patible with equations (6) and (7). The metric is then 

A 2 ( 2 ®\ 2, ,2 , dr 2 2 , 2 

As = - exp T c di + — — + r d 12 

\c 2 J 1 - k M/r 

dfl 2 ee AO 2 + sm 2 6d(f> 2 , (12) 

where $ = $(r), M = M(r) and k = 8nG/c 2 , so that M 
is given in mass units. The source of equation (12) is the 
perfect fluid momentum-energy tensor 

T ab = T: b B + Tq = (e + A) u a u b + (p- A) h ab , 



u a = exp(-$/ C 2 )<5%, 



(13) 



with e, p given by equations (2)-(4). Einstein's field equa- 
tions together with the balance law T^k ; b — for equation 
(12), equations (13) and (17) now take the form 



M' = ( £ + A)r 2 /c 2 , 



k M + (p- A)r 3 /c 2 
2 r (r - itlj : 
$' 

-(e + P)— , 



(14) 
(15) 

(16) 



where a prime denotes derivative with respect to r. Since 
the 4-acceleration is given by u a = <§>'!<? 5„, the constraint 
(7) takes the form 



C 2 ~ P 



= p c exp 



(17) 



where p c = P(0) and $ c = $(0). The subscript c will 
denote henceforth evaluation along the symmetry center r — 
0. Inserting equations (2), (3) and (4) into equation (16) and 
using condition (17) yields after some algebra 



= Co = const., 



(18) 



!L = ^ IJL 

/i P mc 2 

so that n, e and p become functions depending only on P (or 
on $ via equation (17)). Using equation (18) and re-scaling 
n with respect to the central rest-mass density p c = mn c , 
the matter-energy density and pressure for the equilibrium 
MB gas become 



2 / i , 3 



p 



pc 



where the rest-mass density is given by 



P_ 

Pc 



a 1 



(19) 



(20) 



However, it is more useful to express these state variables 
and the field equations themselves in terms of the central 
velocity dispersion (e.g. Binney & Tremaine 1987) 



c 



and the dimensionless variable 



T 2 ' 



(21) 



(22) 



which in the Newtonian limit becomes the "normalized po- 
tential". Combining equations (17), (21) and (22), we obtain 
the "redshifted" temperature provided by Tolman's law 

T = T c exp (^=^) = T c exp *\ , (23) 
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showing the specific temperature gradient of an MB gas in 
thermodynamic equilibrium. 

With the help of equations (17) and (20), we obtain 



— = exp 

Pc 



|^«+a('-- < ~' 4 "*) 



so that e and p become 



= pc 



3a 



2 - 

pa c e 



(<T c /c) 2 * 



(24) 

(25) 
(26) 



The field equation that we have to solve are then equations 
(14) and (15), with s, p and A given by equations (10), (25) 
and (26). However, it is convenient to express these equa- 
tions in terms of equations (21), (22) and the following di- 
mensionless variables 

4nG p c 



M 

leading to 

dM 

dx 

d* 

da; 



— . r o 
ro 

M 



(27) 



3cr 



1+ 2^ 6 



+ A 



M + 



(p/ P c)(a 2 /c 2 )e^^ - A 



x [x-2(a 2 /c 2 )M] 



where p/p c follows from equation (24) and 
A 



q a P^ = 2.538 x lO-^A^Me/EL 

Pc pc 



(28) 



(29) 



while the length scale ro is one third of the "King radius" of 
Binney & Tremaine (1987) and can be given numerically as 



ro 



4.2303 x 10" 



km/sec 



M Q /pc s 

pc 



1/2 



Mpc. 



(30) 



Now equation (28) is a self-consistent system of differen- 
tial equations that can be integrated once numerical values 
for the independent parameters (p c , a c , Oa, h) have been 
supplied, while regularity conditions at the center require 
Mc = M'c = * c = * c = 0. Notice that by setting Q. A = 0, 
we obtain the general relativistic equations of hydrostatic 
equilibrium for a Maxwell-Boltzmann gas characterized by 
the equation of state of equations (2) and (3)— (4). 

An alternative system to that of equation (28), that is 
easier to handle numerically, follows by expressing p, e, p as 
functions of (3 by means of equations (19) and (20), so that 
(3 replaces * in equation (28). However, the variables M 
and in equation (28) are more convenient to explore the 
Newtonian limit. 



5 THE NEWTONIAN LIMIT VS. 

NEWTONIAN CONFIGURATIONS 

We can identify Newtonian conditions, within the framework 
of General Relativity, as the weak field limit of the theory 



for the metric of equation (12) and source of equation (13). 
By looking at equations (23), (24), (25)- (26) and (28), di- 
mensionless variables \& and M (which are not, necessarily, 
small) always appear multiplied by the ratio a 2 /c 2 . Since a c 
is a characteristic velocity of the MB gas, we can examine 
Newtonian conditions by expanding all relevant quantities 
in powers of this ratio. Using relation (27) we find to first 
order in a 2 / c 2 



exp 



1 



2 (j> - $ e 

c 2 

kM 



c 1 x 



(31) 



Thus, the Newtonian limit of General Relativity for the 
sources under consideration is characterized by 

°- « 1 =► ± « 1, ^ « 1, (32) 
c c 2 r 

so that gtt ~ —1, g rr ~ 1 and we can identify (in this limit) 
the Newtonian potential with $. Also, at first order in o 2 /c 2 
we have 

.2 „2 

1 -C 



1 21 * + 

2 c 2 a 2 



-( CTc / c ) 2 * 



-(<Tc/c) 2 * 



and so equations (23), (24) and (25)-(26) become 

2 



P_ 

Pc 



pc 



1 - 



= T c 



■ £ < e, 



In 



(33) 



(34) 



thus, a spherical MB gas in its Newtonian limit approaches 
the "isothermal sphere" with T = T c (or, equivalently a = 
(T c ) for all x (see figure 1). The field equations (28) become 
(at first order in a c /c ) the system 



dM 

dx 


[••( 


-IS 


d* 


M 




da; 


X 2 





+ A 



- A 



x, 



(35) 



where we have kept, together with rest-mass energy den- 
sity (oc p), the other two contributions to the total matter- 
energy density: internal energy ( oc p oc a 2 p) and that of the 
A field (A), as we do not know a priori which one of the latter 
two contributions is larger. In the case without the A field 
(A = 0), we only need to compare rest-mass and internal 
energy densities. From equations (32) and (33), it is evident 
that in this case the Newtonian limit of equation (35) can 
only lead to the equilibrium equations for the Newtonian 
"isothermal sphere" 

dM 
dx 
d£ 

da; x^ 
whose corresponding Poisson equation is 
2 d*\ * 



* 2 

— ex, 



M 



(36) 



-— — (a 

x 2 dx \ 



dx J 



(37) 
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However, if we have the contribution of the A field (besides 
the internal-energy term) then, before deciding which terms 
can be neglected, we should verify first if A is comparable 
in magnitude to the internal energy density term that is 
first order in ct 2 /c 2 . This comparison is especially important 
in the "force" equation (i.e. the second of equations (35)) 
where this term and A have opposite signs. 

For a wide variety of galactic objects (from dwarf 
galaxies to clusters), we have the range 4km/sec < a c < 
1500 km/sec (e.g. Kleyna 2002, Wilkinson 2002, Salucci & 
Burkert 2000), hence 

2 

10~ 9 < % < 2.5 x 10~ 5 , (38) 
c 2 

On the other hand, for currently accepted estimates of cen- 
tral halo density: 0.001 M /pc 3 < p c < 1 M /pc 3 , (e.g. 
Firmani et al. 2001, Shapiro & Iliev 2002, Dalcanton & 
Hogan 2001) we have 

(39) 



4xlO~ 6 fi A h 2 < A < 4xl0~ 3 Ov/i 2 . 



Therefore, since J1a — 0.7 and h ~ 0.65, it is evident that 
near the center we have p/p c = e* w 1 and so ct 2 /c 2 and 
A can have comparable magnitude (especially for structures 
with large p c ~ IMq/pc 3 ). However, as x increases the 
A term remains constant while p/p cr it decreases, implying 
that for large x the constant A term will always end up 
dominating. Also, since the Newtonian potential is negative 
with |$ c | > |$| so that $ < 0, we can safely assume that 
A 3> (ctc/c) 2 e* holds even for small values of x. 

Therefore, we can ignore the contribution from hydro- 
static pressure in equations (35) and replace them by 

AM 

Ax 
d* 
da; 
so that 



(e* + A) 



= T + Ax, 

X 1 



(40) 



_L A. 

x 2 Ax 



(. 



d* 



da; / 



2 A. 



(41) 



These equations are the Newtonian limit of equation (28), 
based on removing only those terms that can be neglected 
for all x in terms of expansions around o c /c and of com- 
parison with rest-mass density. However, it is interesting 
to remark that if we start from the Newtonian isothermal 
sphere, as given by equations (36) and (36), and within a 
purely Newtonian framework (ie not the Newtonian limit of 
General Relativity) we generalize it by adding a A-type re- 
pulsive interaction, we would not obtain equations (40) and 
(41), but the set 



AM 

Ax 
d* 

da; 



* 2 
= e x , 

M J- x 
-r- = 7 + A £■ 



l_ d_/ 2 d*\ = e *_ 3A] 

x 2 Ax V da; / 



(42) 



(43) 



since in Newtonian gravity the A field only enters into the 
"force" equation in (42) but does not enter in the evaluation 
of the gravitational mass M that follows from integrating 
equation (42). On the other hand, in the Newtonian limit 



equations (40) the form of M (and thus the forms of and 
all variables related to M and ^S) is affected by the presence 
of the A term. 

This difference between a purely Newtonian treatment 
and the Newtonian limit from General Relativity is a signif- 
icant feature certainly worth remarking, since it concerns M 
which is the effective gravitational mass affecting the evo- 
lution of test observers and light rays. Since General Rela- 
tivity is currently the best available theory of gravitation, 
we suggest that the system of equations (40)-(41) should be 
preferred (at least conceptually) to the purely Newtonian 
system of equations (42)-(43). In fact, these systems have 
different solutions and if the differences in these solutions 
are not negligible, this can (in principle) be relevant in dark 
matter tracing and gravitational lensing methods used for 
estimating the properties of dark matter haloes of galactic 
clusters. The long numerical integrations of cosmological N- 
body simulations imply that even minor correction factors, 
of the type discussed here, could have important effects on 
the final result. Figures 2a and 2b illustrate how the correct 
Newtonian limit of a Maxwell-Boltzmann gas in equilibrium 
with a A-field is furnished by equations (40)-(41), just as the 
Newtonian limit without A is the isothermal sphere. How- 
ever, before solving equations (40) and (42) numerically it 
is not possible to know if the differences in their solutions 
are significant or negligible. We will examine this issue in 
section IX (see figures 2 and 3). For the remaining of the 
discussion we will consider only the Newtonian limit (NL) 
equations (40). 



6 EXISTENCE AND STABILITY OF 
CIRCULAR GEODESIC ORBITS 

The relativistic generalization of the "rotation velocity" 
along circular orbits is the rotation velocity of test observers 
along circular geodesies in the spacetime, equation (12). 
Bounded timelike geodesies of equation (12) are character- 
ized by (Chandrasekhar 1983) 8 = tt/2 and by two con- 
stants of motion, "energy" Eq = 2e 2 *^ c i and "angular 
momentum" Lo = r 2 ip, where the dot denotes derivative 
with respect to proper time of the observers. Radial motion 
is governed by 



a 2$/c 



1 - kM/t 
where 



+ V(r) 



V(r) 



2*/c 



is the "effective potential" . The conditions for the existence 
of circular geodesic orbits with radius r = r* are r = r = 0, 
leading to 



a 2$/c 2 



1 — r 9' jc 2 
so that 

$' > 0, l-r$'/c 2 > 0, 



1 $'/c 2 



l-r$'/c 2 ' 



(44) 



(45) 



Stability of these orbits requires r = r* to be a minimum of 
the effective potential. Hence we must have V'(r*) = and 
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V"(r*) > 0, leading (with the help of equations (44)) to 

3 + r $" - 2r $' 2 /c 2 > 0. (46) 

Using the full relativistic equilibrium equations (28) and ex- 
pressing conditions (45) in terms of the variables introduced 
in section IV, the conditions for the existence of circular 
geodesic orbits are, up to first order in cr 2 /c 2 , 

_3 



A x 



M 



xe H 2 XM 



% < 0, 
c 2 



1 + • 



A a; 3 



M a_ 

6 



> 0, 



(47) 
(48) 



while the condition (46) for stability becomes 



(3A 



M 



r 

c 2 I 



4A 2 i 3 + (i 3 e -2M)A 



+ 



(* 2 - 3* - 11) x - 6M 



-3M 



* 2M 

x 3 



<0, 



(49) 



In the Newtonian limit we can neglect the terms multiplying 
cr 2 /c 2 , hence condition (48) becomes trivial, while conditions 
(47) and (49) respectively become the Newtonian limit con- 
ditions for existence and stability of circular geodesic orbits 



M - Xx > 0, 
M - (3A-e*)a; 3 > 0. 



(50) 
(51) 



These conditions coincide with those that would have re- 
sulted had we used the Newtonian limit equations (40) (in- 
stead of the relativistic equations) in equations (45) and 
(46). If we had used the Newtonian equations (42) we would 
have obtained the same existence condition (50) but the sta- 
bility condition would be condition (51), modified by having 
the A term multiplied by a factor of 4, instead of 3. Figures 
3a and 3b depict graphically conditions (50) and (51) from 
the numerical solution of conditions (40) and (42). 

The actual velocity of observers along circular geodesies 
follows from evaluating d<p/dt = <p/t, which with the help 
of the previous equations (see Cabral-Rosetti et al. 2002) 
yields 

2 d* 
ax 



M 2 
Xx 

x 



(52) 



where we have used equation (40) to eliminate d* /dx. Since 
$ becomes the Newtonian gravitational potential in the 
Newtonian limit, circular geodesies orbits exist and v(r) is 
well defined as long as the gravitational force field is attrac- 
tive (ie $' > or equivalently d^/dx < 0). 

For physically reasonable spherical configurations, we 
must demand that existence and stability conditions (50) 
and (51) (or (59) and (60)), as well as that v 2 > 0, should 
hold at least in a finite range < x < x\ around the 
symmetry center, where x \ is a radius where one of these 
conditions breaks down. A sufficient condition for this fol- 
lows by looking at the behavior of conditions (50) and (51) 
around the symmetry center. Expanding M around x = 
with the help of equation (40) and the regularity condi- 
tions M(0) = M x (0) = *(0) = # x (0) = 0, where the 
subindex x denotes derivative with respect to x. We obtain 
M w {l/6)M xxx (0)x 3 = (1 + A):r 3 /3, hence for x « 
both conditions (50) and (51) yield (1 - 2A);r 3 /3 > 0, so 



that the central density must satisfy 



p c > 2 Q\ pc 



(53) 



Had we used the Newtonian equations (42) instead of equa- 
tion (40), we would have obtained p c > 3f2Ap C rit, repre- 
senting a very small departure from the Newtonian limit 
expression (53). 

The fact that physically reasonable configurations re- 
quire a minimal bound on p c is an interesting result that 
has no equivalent in the isothermal sphere and the King ha- 
los. Since e* = p/p c and setting A = in equations (50), 
(51) and (52), it is evident that stable geodesic circular or- 
bits, of any radius and for any value of p c > 0, always exist 
for the isothermal sphere and the King models (see section 
VIII) as long as M and p are positive. If there is a A-field 
then for very small central densities violating condition (53), 
conditions (50) and (51) do not hold in a range around the 
symmetry center, though these conditions do hold for a set 
of ranges like < Xi < x < Xj (see figure 3c). However, 
these cases with p c < 2 pcrit are unphysical and will not 
be examined any further. 

Bearing in mind equations (10) and (29), condition (53) 
is satisfied by all known virialized galactic structures whose 
current estimates of p c are between three and six orders of 
magnitude larger than p crit . In fact, the extremely low values 
p c < 2Qa Pcrit would apply today only to large superclusters, 
still "in-falling" within a linear regime, whose density con- 
trast is of order unity. Obviously such large scale structures 
cannot be equilibrium configurations (hence the remarks in 
Nowakowski 2001, Nowakowski et al. 2002 are obvious for 
present day structures). Moreover, p crit is a (model depen- 
dent) function of the cosmic era, hence the minimal p c for 
stable configurations varies with the virialization redshift 
z viT . Assuming a A-CDM model for which Qcdm + = 1, 
the critical density at any z is given in terms of p„ it today 
by 

p„ lt {z) = [(1-O0(1 + z) 3 + Oy] Pcnt(z = 0). (54) 

Thus, for dwarf galaxies which virialized at large z (i.e. 
Zvir < 30), equations (53) and (54) imply 



p c (z = 30) > 1.8 x 10 4 p crit ( 2 = 0) ~ 0.008 M /pc 3 , 



(55) 



a minimal bound on p c that agrees with current estimates 
(Firmani et al. 2001, Shapiro & Iliev 2002, Kleyna et al. 
2002) for central halo densities of these galaxies. Within the 
current hierarchical models of structure formation, where 
large structures form as the result of the merger of smaller 
ones, given Liouville's theorem, it is clear that the cen- 
tral densities of the first bound structures will also repre- 
sent an indicative upper bound on the central densities of 
larger structures. It is interesting that dynamical studies 
across wide ranges of galactic haloes consistently find cen- 
tral densities above the stability requirement at z w 30 of 
0.008M Q pc~ 2 . Figure 4 illustrates the minimal bound on p c 
given by condition (55). 
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7 CUT OFF SCALE OF A AND COMPARISON 
WITH OTHER LENGTH SCALES 

Assuming that condition (53) holds, each one of conditions 
(50) and (51) will necessary break down at some radii, 
xi, X2, given by the smallest positive roots of the equations 

M(xi)-\xl = 0, (56) 
M(x 2 ) - [3A-e* ((a;2) ]4 = 0. (57) 



These radii provide "cut off" length scales that do not exist 
when A = 0. We will compare these length scales with other 
physical cut off scales, such as the virialization radius and 
the radius for the onset of Antonov instability. 

It is useful to express conditions (50) and (51) in terms 
of the volume average of rest-mass energy density p. Un- 
der Newtonian conditions the total matter-energy density 
is Etotai = (p + ^Apcrit) c 2 , while its volume average is 



(^total) 



4vr J r (p + fiAp C rit)r 2 dr _ 3M 



c 2 4n Jq r 2 dr r 3 

= <p) + Ovp„ lt , (58) 
therefore, conditions (50) and (51) can be given respectively 



<P> 
1 / \ , 3 



> 2^APcrit, 

> 2fiAPcrit, 



(59) 
(60) 



where we used equation (29) and M/r 3 = p c M/x [i . Equa- 
tions (56) and (57) then become 

(p)(a;i) = 2 fU pent, (61) 

\(p){x2) + \p(x 2 ) = 2fi A p crlt , (62) 

From the volume averaging definition (58), we obtain by 
integration by parts 

Hence, for < x < x\, where x\ — roX\ is either one of 
xi, X2, we have have: (p) > p. Since condition (59) holds 
in this range, then p x < follows from ty x < (because 
px oc fy x exp(\E')). Therefore for any x in the range of in- 
terest we have (1/4) (p) + (3/4) p > (p). Hence, the root x 2 
corresponds to larger density values than xi, but since both 
p and (p) decrease with x, we must necessarily have 

X2 < xi, (64) 

which is a reasonable result, since it is to be expected that 
the repulsive A-field will destabilize those circular orbits 
with radii X2 < x < xi that lie near the length scale xi 
beyond which no such orbits exist. As we show in section 
IX, the range X2 < x < x\ corresponds to a transition zone 
between the isothermal region < x < X2 and the asymp- 
totic region x > x\ dominated by the A-field. Thus, since 
x = X2 marks the maximal radius of a stable circular orbit, 
we can consider this radius as a "cut off" scale for isother- 
mal spheres in hydrostatic equilibrium with a A-field. In the 
remaining of this section we will use simple qualitative argu- 
ments to compare other cut off length scales of the isother- 
mal sphere with x\. Since numerical estimates (see figure 
5a) show that X2 is about two thirds of x\ (see section IX) , 
this comparison can be extended to x 2 as well. 



Consider the virialization radius, r vir , the currently ac- 
cepted physical cut off scale associated with the virialization 
of galactic structures (Padmanabhan 1995). The condition 
to determine the value of this length scale follows from the 
Newtonian spherical "in fall" (or "top hat") model as Lokas 
& Hoffman (2001), Iliev & Shapiro (2001) and Padnamab- 
han (1995) show, 



(p)(r vir )| Zvir m Ap crlt (z vir ), 

or, equivalently (using equations (27) and (58)), as 



M(x vir )\ 2 



A +n A Pcrit.(Zvi,) 
3 Pc(Zvir) 



(65) 



(66) 



where the factor A is a model dependent (see Lokas & Hoff- 
man 2001 for a discussion on this issue) contrast factor be- 
tween the average rest-mass density of the overdense region 
and the critical density of the cosmological background, both 
evaluated at a virialization epoch, < z viT < 30, which 
depends on the galactic structure and/or a specific struc- 
ture formation model. From equations (61) and (65), we can 
see that both (p)(ri) and (p)(r vir ) are proportional to p„ it . 
Therefore, for all z < z vir , once the galactic structure has 
virialized and the equilibrium equations are valid, we have 



(P)(ri) 
<P>(rv it 



2Q A 



(67) 



or the same expression with a factor 3 instead of 2 if we had 
used the Newtonian equations (42). Since it is reasonable to 
assume (to a good approximation) that (p) ~ p ~ r~ 2 for 
r < ri, a rough qualitative estimate of the ratio ri/r vir is 



ri 

fvir 



2fi A 



1/2 



(68) 



Considering that A ~ 100 in the case when the cosmological 
background complies with fi CDM +fiA = 1 (Lokas & Hoffman 
2001), this estimate means that for present day virialized 
structures ri is approximately one order of magnitude larger 
than r vir . Because of equation (64), r 2 should be less than an 
order of magnitude larger than r vir . The precise numerical 
value of x vir is displayed in figure 5b, while the qualitative 
estimates mentioned before are corroborated numerically in 
figures 6 and 7, showing that T2 is about seven times larger 
than r vir . Using the Newtonian equations (42) yields smaller 
n and r 2 by a factor of y/2/3. In either case, these results 
indicate that the A-field should have a negligible effect on 
the virialization process. 

Another physically motivated length scale, character- 
istic of the isothermal sphere, is the critical radius asso- 
ciated with the Antonov instability, or with the onset of 
the gravothermal catastrophe. This radius follows from de- 
manding that the structural parameters of an isothermal 
sphere should correspond to a stable thermal equilibrium, 
characterized by a second variation of a suitably defined 
entropy functional. As discussed in Padmanabhan (1990), 
for all isothermal spheres extending up to a critical ra- 
dius R so that E is the Newtonian total energy at R, 
the constraint RE/{GM 2 ) > -0.335 must hold, though 
stable configuration with local maxima must also satisfy 
p c < 709 p(R), while isothermal spheres complying with 
RE/(G M 2 ) > -0.335 but for which p c /p(R) > 709 cor- 
respond to metastable configurations in which the entropy 
extrema is a saddle point. The critical radius R complies 
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with 

Pc 

p(R) 



709, 



(69) 



corresponding to the maximal stable isothermal sphere for 
which entropy is a local maximum. Using equation (29), wc 
can find the ratio 



p(R) 



5.55 x 10 3 
h 2 



Pc 



(70) 



Pcrit ft" M /pC 3 ' 

which can be compared to the ratio of (p)(r vir ) to p crit given 
by equation (67). This yields approximately 



(P)(r 



P(r„ 



h 2 A M /pc 3 



p(R) 



p(R) 



5.55 x 10 3 



(71) 



but, if R < n, then we should have approximately (p) 
p ~ r~ 2 for r < R, hence we obtain a ratio 

R hA 1 / 2 /M /pc 3Xl/2 



74.16 



Pc 



and also the following ratio analogous to (68) 

1/2 



n 
R 



74.16 



(2 0a) 1 / 2 ft V M ©/P c3 



(72) 



(73) 



Assuming that h — 0.65, Qa = 0.7, equations (72) and 
(73) show that for values (consistent with current estimates) 



around p c 



10" 



10 M /pc (see figures 5 and 6), the 



virialization radius is about the same order of magnitude 
as the length scale R of equation (69) and about one to 
two orders of magnitude smaller than ri . An estimate of the 
value of R in Mpc follows from the fact that the variables 
(27) that we use coincide (save notation differences, a minus 
sign in \1/) with those defined in equations (4.21) and (4.22) 
of Padmanabhan (1990) which yield x{R) = R/r = 34.2. 
Using equation (30) we obtain 



r(R) 



1.467 x 10" 



km/sec 



M /pc ; 

Pc 



1/2 



Mpc. 



By looking at the numerical solutions for the equilibrium 
equations (40) and (42), we will verify in section IX the 
qualitative estimates of ri,r2,r vir and R that we have pre- 
sented above. These estimates tend to support the conclu- 
sion that the effects of a A field can be ignored in the study 
of virialized structures. 



8 THE KING MODELS 

It is interesting to compare the effect of a A-field on a MB 
distribution with other distribution functions (but without 
A). We examine briefly the case of the Michie-King models, 
obtained by truncation of the MB distribution at a given 
maximal escape velocity (Binney & Tremaine 1987). This 
distribution functions leads to equilibrium equations similar 
to (36), but with p/p c given by Binney & Tremaine (1987) 
and Katz (1980) 



P_ 

Pc 



= * + *o, 



*0 



(75) 
(76) 



-i log 10 x 
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Figure 1. Each pair of curves provides the radial temperature 
variations within a MB sphere characterized by a given central ve- 
locity dispersion (values indicated in km/sec). The curves showing 
monotonous decay correspond to A = 0. In the case with nonzero 
cosmological A term, the temperature curves stop decreasing and 
oscilate around a constant temperature (see text). It is evident 
that as the velocity dispersion of the halo decreases, the halo 
tends to an isothermal configuration, valid to a high degree in 
any astrophysical structure, <r c < 1000 km/sec. 



K{X) 



e x erf(A 



l/2x 



2X 1 / 2 



(77) 



where ^ is given by equation (22), while $6 = &(xt) for 
x = x b satisfying the constraint K(xb) = => p(x b ) = 0. 
The Xb is known as the "tidal radius" marking a theoretical 
cut off scale whose value depends on ^o- The equilibrium 
equations for the King models expressed in terms of the 
dimensionless variables x and M. are then 



(74) 2^1 



AM 

dx 
d* 

da; 



2 Km 

M 



(78) 



Since *(0) = V&o and M(0) = 0, the integration of these 
equations requires one to specify a given value of ^o, leading 
to a family of numerical solutions parametrized by \&o- The- 
oretical studies of the King distribution applied to globular 
clusters (see Katz 1980) yield the stability range $o < 8, 
though there are observations and theoretical estimates al- 
lowing for values up to *o ~ H- 

As we will show in the following section, the dynami- 
cal variables of the King models behave asymptotically in 
a very different manner from those of the MB models with 
A. However, the cut off scale r\ in the latter is of the same 
order of magnitude as the tidal radius of King models with 
deep potential wells (^o ~ 10). 



9 NUMERICAL ANALYSIS 

We solve numerically the systems of equations (28), (36), 
(40), (42) and (78), assuming the regularity conditions 
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(a) 




Figure 2. a) Relative potential vs. radius, for MB spheres in the 
absence of a A term. The curves are characterized by sequentially 
decreasing central velocity dispersions, ranging from relativistic 
values (upper curves) to newtonian values (lower curves). The 
solid curve at the bottom corresponds to the classical isother- 
mal sphere solution. As the central velocity dispersion decreases 
the sequence accurately converges to the isothermal sphere, b) 
Same as panel (a), but in the presence of a cosmological A term. 
The curves converge to the solid curve marked NL associated with 
equations (40), hence it is clear that this is the correct Newtonian 
limit for the MB spheres obtained from the full General Relativity 
treatment. The curve marked N corresponds to the naive Newto- 
nian limit associated with equation (40). It is evident that this is 
not the correct newtonian limit of the MB sphere when there is 
a A-field. 

M c = M' c = * c = *' c = 0, for Q A = 0.7, h = 0.65 and 
suitable values of the free parameters <t c , pc, ^o- 

9.1 The isothermal limit with a A field 

It is interesting to examine the Newtonian limit of the MB 
gas (with and without the A-field) by looking at how an 
isothermal behavior emerges as temperature gradients de- 
crease for sequentially decreasing values of the ratio cr c /c, 
starting from relativistic values down to smaller Newtonian 
velocities. Using from the numerical solution of equa- 
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Figure 3. Values of the existence (M — Ax 3 ) and stability 
(M — (3A — e*)x 3 ) conditions for circular orbits vs. the dimension- 
less radial coordinate x for MB spheres with p c = 0.01Mq/pc 3 . 
Panels (a) and (b) respectively depict the cases of the Newtonian 
limit and the naive Newtonian treatment. In the latter case the 
existence and stability radii are about two thirds smaller than 
those of (a). Panel (c) For depicts the case where central density 
is lower than the threshold p c > 2f2>,. In this case no central 
structure can exist, as the existance and stability conditions are 
only met in disconnected ranges not including the origin. 
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tion (28), we plot in figure 1 the normalized temperature, 
T/T c , given by equation (23) vs log 10 x for four values in 
the range a c /c — 0.08 down to a c /c — 1.6 x 10~ 2 , that 
is from a c — 25000 km/sec to a c = 5000 km/sec (still a 
dispersion velocity much larger than that of any galactic 
structure). As shown in this figure, the temperature gradi- 
ents of the MB gas are already small (T/T c < 0.2) even 
for a c = 25000 km/sec, suggesting that even for more rela- 
tivistic values of a c /c the linearized expression (34) can be 
used instead of equation (23). These gradients steadily de- 
crease for smaller u c /c, becoming negligible for the velocities 
characteristic of galactic structures, and so approaching in 
the limit a c /c — > an isothermal condition T = T c . In the 
case A = the temperature curves steadily decrease so that 
T/T c — » as x — * oo. However, if there is a A-field, then 
these curves stop decreasing at x — x\ and begin to oscillate 
around a constant value, T = T\, that can be identified as 
a characteristic temperature of the A-field. 



9.2 The correct Newtonian limit 

The Newtonian limit of equation (28) can also be illustrated 
by means of sequences of curves corresponding to $ obtained 
from equation (28) for sequentially decreasing values of cr c / c. 
The curves of figure 2a (A = 0) clearly tend to the thick curve 
at the bottom, marked as "IS", corresponding to ^ for the 
isothermal sphere, obtained from equation (36). In figure 
2b we consider the case with nonzero A given by equation 
(29). It is evident that the sequence of curves tends to the 
thick curve, marked NL, corresponding to it obtained from 
equation (40) and not to the thick curve, marked "N" , that 
corresponds to ^ obtained from equation (42). Hence, the 
correct Newtonian limit of a MB gas with a A field is given 
by equations (40)-(41) and not by the Newtonian equations 
(42)-(43). 

9.3 The "cut off" length scales 

Numerical solution of the systems (40) and (42) allow us 
to depict graphically conditions (50)-(51) (figures 3a and 
3b), conditions (59)-(60) giving the same results, where the 
length scales for the existence and stability of circular or- 
bits, x\ and X2, follow from the roots given by equations 
(56)-(57) and (61)-(62). Assuming p c = O.OlM /pc 3 , the 
following numerical values emerge from figure 3a: si fs 800, 
X2 ~ 555, while for the Newtonian equations (42) we obtain 
smaller values (see figure 3b): x\ w 660, X2 ~ 460. For the 
same value of p c and for the Newtonian limit equations (40) 
we obtain from figure 5: x vir w 72 and x(R) m 35. These 
values of x can be translated into actual distances for spe- 
cific galactic structures by using r = ro x where ro is given 
by equation (30). The resulting length scales are shown in 
Table 1. 

As we mentioned in section VI, the conditions (50)- 
(51) for the existence and stability of circular geodesic or- 
bits are violated for p c < 2 Qa pent ■ This possibility is il- 
lustrated by figure 3c, where these conditions are plotted 
for p c — £7a p cr it, clearly showing how conditions (50)-(51) 
do not hold around the center (though they do hold in dis- 
connected ranges of x that exclude the center). Since both 
Pc and p crit are functions of the cosmic era, condition (53) 
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Figure 4. The solid curve shows the minimum central density 
for stability as a function of virialization redshift. The value at 
z = 30, corresponding to p c = O.OIMqpc -2 , is of the order of 
observed values in galactic structures which can be thought of as 
having formed at around that time or less. 
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Figure 5. Logarithmic plots of the density and average density 
vs. dimensionless radius x, both normalized by p C rit- The inter- 
sects with the horizontal dotted lines define the onset of Antonov 
instability (p = p c /709) and the virial radius (< p >= Ap cr it). 

with the help of equation (54) yields a minimal bound of 
Pc for virialized structures complying with conditions (50)- 
(51). This is illustrated in figure 4, showing that struc- 
tures of relatively recent virialization easily meet this mini- 
mal bound, though for galaxies which virialized at z ~ 30, 
the fulfillment of condition (53) is non-trivial, leading to 
p c > 0.008 Mq/pc 3 , a value that is consistent with current 
estimates of p c . 

In figures 3a-b and 5 we considered p c = 0.01 Mq/pc 3 . 
In order to find out the dependence of the values of xi, X2 
and R on p c , we plot in figure 6 the ratios of these length 
scales with respect to x vil , clearly emerging from this figure 
that xi and X2 have a very weak dependence on p c and are, 
respectively, about 10 and 7 times larger than a; vi r (for the 
Newtonian case these figures are about two thirds smaller). 
However, R is quite sensitive to p c , with x viT ~ R for densi- 
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Length scale 


Dwarf 
galaxy 


Large 
galaxy 


Cluster 


Existence of circular 
orbits, n 


140 kpc 


5 Mpc 


32 Mpc 


Largest stable circular 
orbit, T2 


90 kpc 


3 Mpc 


22 Mpc 


Virial radius, r v i r 


14 kpc 


300 kpc 


3 Mpc 


Antonov instability, R 
p c = 0.01M Q /pc 3 


8 kpc 


140 kpc 


1.3 Mpc 


Antonov instability, R 
p c = 0.001 M Q /pc 3 


15 kpc 


330 kpc 


3.2 Mpc 



Table 1. Length scales for galactic structures. Dwarf and large 
galaxies and clusters are respectively characterized by cr c = 
3, 200, 1000 km/sec. The length scales n, r<i and r vir were eval- 
uated for p c = 0.01 Mq/pc 3 . See also figures 6 and 7. 

ties in the range 0.001 — 0.01 Mo/pc 3 . This is a nice result, 
illustrating that these two physically motivated length scales 
approximately coincide for values of p c close to current esti- 
mates. All length scales we have considered show a clear lin- 
ear dependence on <r c (from equation (30)). This is shown in 
figure 7, displaying the graph of each length scale as a func- 
tion of (7c for a wide range of values of p c (10~ 5 — 1 M Q /pc 3 ). 
Since the dependence on p c is very weak, the strip of lines 
for all values of p c simply makes the lines thicker. 

9.4 Density profiles and rotation curves 

The effects of the A-field can be clearly illustrated by plot- 
ting relevant physical variables of the isothermal halos, such 
as the density profiles and rotation velocity along circular 
orbits. 

Figure 8 depicts a logarithmic plot of the rest-mass den- 
sity, p, normalized by p crit for p c — 0.01 M@/pc 3 , in compar- 
ison with the rest-mass density of the isothermal sphere and 
of King halos associated with various values of (only the 
Newtonian limit "NL" case is shown). It is evident that in 
the range < x < X2 where stable circular orbits exist, 
the density curve of the NL case is almost identical with 
that of the isothermal sphere (marked as "IS"). However, 
in the asymptotic range x > Xi, these curves are com- 
pletely different, since for the IS p simply decays as x~ 2 
while in the NL case it oscillates around a value close to 
p crit . Since p — pee 9 , these oscillations mark relative max- 
ima and minima of the Newtonian potential, hence x = x\ 
marks an unstable equilibrium where this potential becomes 
repulsive (though furhter out it alternates from repulsive to 
attractive). Clearly, this behavior of the asymptotic region 
is dominated by the A-field and is in stark contrast with the 
behavior of the isothermal region < x < X2- The region 
X2 < x < Xi, where circular orbits exist but are unstable, 
is a transition zone between the isothermal and asymptotic 
zones. However, in the physical range < x < x vir , the ef- 
fects of the A-field are practically negligible. Notice that p 
for the NL case is also very different from the density curves 
of the King halos, though the cut off boundary (the tidal 
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Figure 6. Dependence of circular orbits existence and stability 
radii, r\ and r2, and onset of Antonov instability, R, all in units of 
the virial radius, as functions of the logarithm of central density. 
It is evident that both rx and r2 depend only very weakly on p c 
and are always well beyond the virial radius. The radius R is much 
more sensitive to central density, being rather small for larger p c . 
Since a heat conduction mechanism (interaction cross-section) is 
required for this mechanism to operate, it probably does not play 
a significant role. 
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Figure 7. The plot depicts the relation between the indicated 
radii vs. velocity dispersion, where we have marked with shaded 
regions the zones occupied by various astrophysical structures. 
The solid lines represent existence, stability and virial radii for 
a wide range of central densities. The dotted lines correspond to 
radii for the onset of Antonov instability for the two indicated 
values of the central density. 

radius) of some of these halos (^o > 10) can be of the same 
order of magnitude (or even larger) than X2- 

Rotation velocities normalized by a c are given by equa- 
tion (52). These "rotation curves" are plotted in figure 9 
for the NL case for p c — 0.01Mo/pc 3 , in comparison with 
rotation curves of the isothermal sphere and of King halos 
associated with various values of ^o- The figure also displays 
the cut off length scales Xi, X2, separating the isothermal, 
asymptotic and transition zones, as well as the physical ra- 
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Figure 8. Density profiles normalized by p crit as functions of 
the dimensionless radius x. Curves labeled by K8 to Kll give 
classical King profiles of the indicated shape parameter, IS is the 
isothermal sphere, and NL corresponds to the newtonian limit 
of an MB gas in the presence of a cosmological A term. Vertical 
dotted lines indicate the virial radius and the radii for stability 
and existence of circular orbits. 
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Figure 9. Rotation velocity curves for the cases shown in figure 
8. The graphs illustrate that within the Isothermal region, NL 
follows IS to within a factor of 1.4, however, within the Transition 
region, the NL rotation curve falls to zero, (see text) 

dius a; V ir. Since this is not a logarithmic plot, the effects of 
the A-field are more noticeable. For example, in the NL case 
the rotation curve is not flat and decreases significantly al- 
ready in the isothermal region, though it is practically iden- 
tical to the rotation curve of the isothermal sphere in the 
range < x < x viT where all visible matter would be lo- 
cated. The rotation curve of the NL case plunges to zero 
for x = Xi, so that for a range of values x > x\ circular 
geodesic orbits do not exist. This is connected to the fact 
that x — xi marks an unstable equilibrium where the New- 
tonian potential becomes repulsive. However, this potential 
oscillates, hence it becomes again attractive and then re- 
pulsive again, all of this in disconnected domains of x for 
x > x\. As a comparison, the rotation curves of King halos 
(evaluated up to their tidal radii) also decrease but never 
reach zero. 



10 CONCLUSION 

We have examined the equilibrium of a MB gas in the pres- 
ence of a A-field constrained by recent observational data. 
Proceeding from a general relativistic framework and care- 
fully taking the Newtonian limit, we obtain the equilibrium 
equations (40)-(41) that are different form those that fol- 
low from a "naive" Newtonian treatment that would simply 
add the A-field to an isothermal sphere (i.e. equations (42)- 
(43)). We have shown that the presence of this repulsive field 
implies that isothermal configurations are now characterized 
by the length scales, x\, X2, defined by conditions (50) and 
(51) and respectively associated with existence and stabil- 
ity of circular geodesic orbits. We have shown that these 
length scales can also be characterized by the constraints 
given by (59) and (60) on the rest-mass density and its 
volume average, hence x% < x\. The fulfillment of the ex- 
istence and stability constraints in a region containing the 
symmetry center x — yields a minimal central rest-mass 
density p c = 2 Qa Pcrit , and thus a minimal central density of 
p c = O.OO8M0/PC 3 for galactic structures having virialized 
at z ~ 30. Since x = X2 marks the radius of the largest stable 
circular orbit, it is reasonable to consider this radius as the 
cut off length scale of isothermal configurations in equilib- 
rium with a A-field. The numerical solution of the equilib- 
rium equations provides precise numerical estimates of all in- 
volved length scales, showing that X2 is about 7 times larger 
than physically motivated cut off length scales of isothermal 
halos, such as the virialization radius x vir and the radius 
for the onset of Antonov instability R. By looking at the 
graphs of the rest-mass density and rotation curves, we can 
identify isothermal (0 < x < X2), asymptotic (x > xi) and 
transition (X2 < x < xi) regions in which these variables 
change from an almost isothermal behavior to an asymp- 
totic regime dominated by the A-field. It is quite noticeable 
that up to the physical cut off scale, x vir , the behavior of the 
MB gas with the A-field is practically indistinguishable from 
that of an isothermal sphere. Therefore, for the purposes of 
studying the equilibrium state of virialized structures, the 
effects of the A-field are negligible and can be safely ignored. 
However, these effects might be important in the study of 
the gravitational clustering of larger structures that have 
already left the linear regime but are not yet virialized, or 
the formation of galactic structures (though the assumption 
of thermodynamical equilibrium might not be applicable in 
these cases). 

Finally, the zero gradient approximation associated 
with the constant A-field is basically a toy model of more 
general quintessence-like sources described by non-trivial 
scalar fields. The equilibrium of a MB gas together with such 
a source could yield interesting features that are absent in 
the A-field approximation. This problem is currently under 
investigation (Niiez, Matos and Sussman, 2003). 
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